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Abstract 

MicroRNA-mediated regulation of gene expression is characterised by some distinctive features that set 
it apart from unregulated and transcription factor-regulated gene expression. Recently, a mathematical 
model has been proposed to describe the dynamics of post-transcriptional regulation by microRNAs. The 
model explains the observations made in single cell experiments quite well. In this paper, we introduce 
some additional features into the model and consider two specific cases. In the first case, a non-cooperative 
positive feedback loop is included in the transcriptional regulation of the target gene expression. In the 
second case, a stochastic version of the original model is considered in which there are random transitions 
between the inactive and active expression states of the gene. In the first case we show that bistability is 
possible in a parameter regime, due to the presence of a non-linear protein decay term in the gene expression 
dynamics. In the second case, we derive the conditions for obtaining stochastic binary gene expression. We 
find that this type of gene expression is more favourable in the case of regulation by microRNAs as compared 
to the case of unregulated gene expression. The theoretical predictions relating to binary gene expression 
are experimentally testable. 

1 Introduction 

MicroRNAs (miRNAs) are a class of small, noncoding RNAs which regulate gene expression in prokaryotes 
and eukaryotes at the post-transcriptional level. They play critical roles in a number of cellular processes, 
such as stress response, developmental transitions, differentiation, apoptosis etc. [TJ[SJ[3]. The mechanisms of 
regulation by small RNAs differ in specific features in prokaryotes and eukaryotes. These are, however, based 
on the common principle of the regulatory RNA base-pairing with the messenger RNA (mRNA) of the target 
gene inhibiting translation and/or promoting mRNA degradation. Various possibilities have been suggested for 
the fate of the inactive complex of the regulatory RNA and the target mRNA once it is formed [H [5J H] : 
(i) the complex has a finite lifetime followed by dissociation into its free components, (ii) the regulatory RNA 
co-degrades with the target mRNA at the same or different rates and (iii) only the target mRNA degrades 
with the regulatory RNA becoming free for further activity [H [7] . A number of functional features of 
RNA- regulated gene expression have been identified so far [4J [6j [7] . One prominent feature is that of a threshold 
linear mode of action, in which the target gene protein synthesis is highly repressed below a threshold level 
of target mRNA production and activated in a linear fashion once the threshold is crossed. Other significant 
characteristics include: suppression of protein fluctuations in the form of translational bursts, rapid response 
times and filtering of transient signals [7] , sharpening of spatial expression patterns [5] and prioritisation of the 
expression of target genes in the case of a single RNA regulating the expression of multiple genes [5]. 

Recently, a mathematical model has been proposed [I] to describe the biochemical interactions and kinetics 
of miRNA-mediated regulation of gene expression. The model has been experimentally validated in single 
cell measurements using quantitative fluorescence microscopy and flow cytometry. The experiments clearly 
demonstrate the high repression of target protein synthesis below a threshold level of target mRNA production, 
and a sensitive response above the threshold. A major finding is that the strength of the repression below 
the threshold has considerable cell-to-cell variation in a population of genetically-identical cells. This indicates 
that stochasticity may have non-trivial consequences in miRNA-mediated regulation of gene expression. Figure 
[T] shows a sketch of the different biochemical processes involved in the miRNA-mediated regulation of gene 
expression, and which form the basis of the mathematical model proposed in pQ. The mRNA, as shown in 
Fig. [TJ is either free (m) or is part of a miRNA-mRNA complex (m*). The target mRNA is transcribed 
with rate constant k m and has a natural degradation rate constant j m . The rate constants k on and k ff 
are associated with the formation and dissociation of the bound complex of mRNA and miRNA. The bound 
miRNA becomes free either by unbinding the target mRNA (rate constant k ff) or by degrading the mRNA 
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Figure 1: A schematic diagram describing the biochemical events involved in miRNA-regulated gene expression. 
The target mRNA is synthesised and degraded with rate constants k m and 7 m respectively. The rate constants 
k on and k j / describe the formation and dissociation of the bound complex m* of free mRNA m and the miRNA. 
The bound miRNA degrades the mRNA with rate constant 7 m . . 



(rate constant j m * ) . The model incorporates the important feature of molecular titration similar to the protein- 
protein titration analysed earlier |10] , Protein sequestration occurs when a repressor protein binds an active 
protein thus forming an inactive complex. As shown in |10) . the regulatory mechanism, termed molecular 
titration, can generate ultrasensitive input-output responses. In the case of miRNA-mediated regulation of 
gene expression, the inactive complex is that of mRNA and miRNA. Molecular titration is responsible for the 
observed sensitive dependence of protein expression on target mRNA input around a threshold level of target 
mRNA production. 

In this paper, we analyse the dynamics of the model, the key biochemical processes of which are shown in 
Fig. [TJ from two different perspectives. In the first case, an additional process representing the autoactivation 
of the target gene expression by the target protein is included. In the second case, a simple stochastic version of 
the model of Fig. Q]is considered. In this model, the only stochasticity is associated with the random transitions 
between the inactive and active expression states of the gene. In the inactive state, there is no transcriptional 
activity, i.e., mRNA synthesis. Transcription is only initiated in the active state of the gene. The degradation 
of the mRNA occurs in both the inactive and active states of the gene. We show that binary gene expression is 
a possible outcome in both of the scenarios described above. In the deterministic case, binary gene expression 
implies bistability, i.e., the coexistence of two stable expression states. In the stochastic case, the distribution 
of mRNA levels is bimodal, i.e., has two prominent peaks. In the following, we analyse the deterministic and 
stochastic models to investigate the origins of binary gene expression. 



2 Deterministic Model 

The total concentration of miRNA is assumed to be constant, since the miRNA turnover is slow compared to 
the timescale of gene expression [4J. The following set of differential equations describe the dynamics of the 
model. 

dm p m p . 

—jr = k m + — — ■ k on m.mi + koffm -7 m m (1) 

dt K + p 

dm* 

= k on m.mi - k ff m - j m , m (2) 



dt 

dp 



fU k p m-7 p p (3) 



The conservation condition for miRNA is 



mix = rni + m* (4) 

In the equations above, mi and mix represent the free and total miRNA concentrations. The second term on the 
r.h.s. of Eq. (1) represents the autoactivation of the target gene expression. The rate constant f3 m represents 
the maximum rate of mRNA synthesis due to autoactivation, with K denoting the equilibrium dissociation 
constant for the binding of the regulatory protein at the promoter region of the target gene. The rate constants 
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Figure 2: (a) Steady state protein concentration p st versus a. The parameter values are kg = 0.01, K = 10, 
<j) — 0.2 and X a = 0.1. (b) Steady state protein concentration p st versus X a . The parameter values are kg — 0.01, 
K = 10, 4> = 0.3 and a = 13. (c) Phase diagram in the <fi — a plane showing a region of bistability between 
two regions of monostability. A continuous path between the two monostable regions is also possible. The 
parameter values are kg = 0.01, K = 10 and A Q = 0.1. (d) Phase diagram in the cf> — A Q plane showing a region 
of bistability between two regions of monostability. The parameter values are kg = 0.01, K — 10 and a = 13. 



k p and 7 P correspond to protein synthesis and degradation, respectively, with p being the protein concentration. 



In the steady state, ^2 = q, = and & = and one obtains the equations 



fc m + ~tt~t~ k on m.mi + k of f m* - j m m = 0, 7* = 7 m » .mi T (5) 



K + p 



m.mi T k off + j m * 
mr = -, \ = — LL (6) 



m + X 



k 

p = —m (7) 

7 P 



From these equations, the steady state protein concentration satisfies the equation 

a p st (j>Pst , a , 

Pst = kg + — — — 8) 

K+p st Pst + X a 

with kg = — , a — — , 4> — \ a = \ — and 7q = ^f^. 

Eq. ([S]) has two stable steady state solutions, i.e., bistability in specific parameter regions. Figure 2(a) shows 
a plot of the steady state protein concentration, p st , as a function of the parameter a. The stable steady 
states are represented by solid lines whereas the dotted line describes the branch of unstable steady states. 
The plot exhibits hysteresis since the discontinuous transitions from the lower to the upper branch, and from 
the upper to the lower branch, occur at different values of the parameter a, the so-called bifurcation points of 
the associated dynamics. Bistability constitutes a universal theme in several cell biological processes |11[ PT2"] 
and signifies that a cell has a choice between two stable expression states for the same parameter values. The 
choice of the specific stable state depends on the previous history of the system. Figure 2(b) exhibits a plot 
of p s t versus the parameter A Q . Figure 2(c) shows a phase diagram in the <j) — a plane in which a region of 
bistability separates two regions of monostability, corresponding to low and high expression states, respectively. 
Figure 2(d) shows a similar phase diagram in the <fi — A Q plane. Bistability, in general, is an outcome of 
dynamics involving positive feedback and sufficient nonlinearity. The latter condition is usually achieved via 
the binding of the regulatory protein molecules at multiple sites of the promoter region of the gene, or when the 
regulatory proteins form multimers, such as dimers and tetramers, which then bind specific regions of the DNA 
(cooperativity in regulation) [T3J [T3]. In the case of transcriptional regulation, Tan et al. [T5] have proposed 
a new mechanism by which a noncooperative positive feedback loop combined with a nonlinear protein decay 
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term are sufficient to generate bistability. This novei type of bistability was demonstrated in the operation of a 
synthetic gene circuit. The circuit contains a single autoregulatory positive feedback loop in which the protein 
product X of a gene promotes its own synthesis in a noncooperative fashion. The protein decay rate is a sum 
of two terms, the natural degradation rate and the dilution rate due to cell growth. In the synthetic circuit, 
the production of X inhibits cell growth so that the dilution rate of X and hence the protein decay rate are 
reduced. This gives rise to an effective positive feedback loop, since an increased synthesis of X proteins leads 
to a greater accumulation of the proteins, which in turn activates further protein synthesis. The combination 
of two positive feedback loops gives rise to bistability in a parameter regime in the absence of cooperativity. 
A related study by Klumpp et al. (16j has also demonstrated the generation of a positive feedback loop due 
to cell growth inhibition by a protein. The nonlinear protein decay term in these cases ('emergent bistability') 
has the same form as in Eq. (8). The origin of the nonlinearity is, however, different in the case of the post- 
transcriptional, i.e., miRN A- mediated regulation of gene expression. The effect of a transcriptional positive 
feedback loop on miRN A- regulated gene expression has been investigated in a number of earlier studies |17j . 
Specifically, bistability has been demonstrated in cases where the transcriptional positive feedback loop involves 
some form of cooperativity. Our present study establishes a new origin of bistability in miRNA-regulated gene 
expression, based on a non-cooperative transcriptional positive feedback loop and a nonlinear protein decay 
term. Since our model closely follows the experimentally-tested model proposed in Ref. [3], the inclusion of a 
transcriptional positive feedback loop in the original gene circuit could provide an experimental test of the new 
mechanism of bistability suggested by us. 



3 Stochastic Model 

The model is a generalisation, incorporating miRNA-mediated regulation of gene expression, of an earlier 
stochastic model |18| in which the effect of stochasticity on unregulated/transcriptional factor-regulated gene 
expression was considered. More specifically, the condition for obtaining stochastic binary gene expression in 
terms of two gene expression parameters, rj = and T2 — ■d L - was derived in the previous study. The rate 
constants k a and kd are the activation and deactivation rate constants for transitions between the inactive and 
active states of the gene. The rate constant "f m is the degradation rate constant of the gene expression product, 
which could be either mRNA or protein. In the present study, we determine the steady state distribution of 
mRNA levels. We assume that the concentration m* of the mRNA-miRNA complex attains its steady state value 
at an earlier time point than the concentration m of free mRNA. As already mentioned, the only stochasticity 
in the model considered here is associated with the random transitions between the inactive and active states 
of the gene. In each state of the gene, the mRNA concentration evolves according to the equation, 

dm 7* m 

— = k m z — - 7 m m (9) 

at m + A 

Eq. © is obtained from Eq. (1) by ignoring the positive feedback term and substituting the steady state 
concentration of m* (Eq. (6)). The variable z = 1 (0) when the gene is in the active (inactive) state and switches 
values with stochastic rate constants k a (0 — > 1) and kd (1 — > 0) . Let pj(m,t) (j = 0, 1) be the probability 
density function when z — j with the total probability density function p(m,t) — po(m,t) + pi(m,t). The rate 
equations for the probability density functions are given by 

dpv{m,t) d r/ 7*m . . . 

57 = --5— {( — r - 7m m)po(m,t)\ + k d pi{m,t) - k a p Q (m,t) (10) 

at am m + A 

dp\(m> : i\ d 7*?Ti 

ttt — = -^—{{k m ; — r - j m m)pi{m,t)} + k a p (m,t) - k d pi(m,t) (11) 

at am m + A 

The first terms in Eqs. (10) and (11) are the so-called 'transport' terms, representing the net flow of the prob- 
ability density. The second terms represent the gain/loss in the probability density due to random transitions 
between the state j (j = 1, 0) and the other accessible state. In the steady state, both and are zero and 
the total probability density function p satisfies the equation 

dp = t^fHp | (/'M-^p 

dm f( m ) 1 — f( m ) 

where 



(12) 



/( m ) = 1 ^fr + Jne w m, (13) 



* 7* 
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The steady state solution for p is given by 
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Figure 3: Plots of p(m) versus m (Eq. (14)) in the steady state [a\ < and (14 + 0,5 < 0). The values of n, 
f2 and A are shown in individual boxes. The other parameter values are ~/ new = 0.1, 7^ el0 = 0.2. (a) n < 1, 
r<z < 1, (b) 7*1 > 1, T2 < 1. For a wide range of A values, binary gene expression is obtained in each case. 



p(m) = Nm ai (d 2 +m) a2 (d 3 + m) a3 (d 5 - m) a4+a5 (d e + m) Q4 " a5 (14) 
where N is the normalisation constant, 

a ^ _ kg . A . ^ ^ 

7neu> 7neu>A 
a ( 7raeu; \ ^ 

a 3 = 2 

fed 1 



a 4 



27 



nci.r 



a 5 = 7— ( ^ ) / ; = ( 15 ) 

k m ilne W 2 ^ Xlnew + ( fe^+^d )2 



and 



d 2 = + A 

1 new 

<fe = A 



o + o V (7^ e «: + ^Inew - 1) + 4A7„ 



e?5 

2 7 

new 27 

(7ne» + A 7«eu> - 1) , 1 



+7y V (inew + X lnew - l) 2 + 4\j new (16) 

Putting 7* eul = and A = 0, i.e., considering only unregulated gene expression, one recovers from Eq. (14) the 
beta distribution QUI [EH]: 



p(m) — N\ m!"<rn _ m)^™ 1 - ) (17) 

7m 

where Yi is the normalisation constant. In this case, binary gene expression is obtained in the parameter 
regime r\ = — and r2 = both < 1. Two prominent peaks in the probability density function appear 
when r\ and are comparable in magnitude. In the case of transcription-factor regulated gene expression, 
the effective activation and deactivation rate constants, k a (s) and k d (s), are functions of the regulatory protein 
(transcription factor) concentration s. The gene expression response to a regulatory stimulus may be either 
graded or binary. The response is quantified in terms of the concentrations of mRN As/proteins. In graded 
response, the average steady state concentration of the gene expression product varies continuously as the 
concentration s of the regulatory molecules is changed, until a saturation level is reached. In the case of binary 
response, gene expression occurs at either of two average levels (say, low or high) and expression at other levels 
is minimal. The fraction of cells in the low/high expression level changes as s is changed. This gives rise to 
a bimodal distribution in the protein/mRNA levels in an ensemble of cells. In the case of transcription- factor 

k ( s) k ( s) 

regulated gene expression, binary gene expression is obtained when n = ° < 1 and r2 = ^ < 1. One 
should point out that in the parameter regime in which binary gene expression occurs, unimodal distributions are 
obtained when r\ << r-i or « r\. The full parameter regime is associated with the system exhibiting binary 
response to changing activation and deactivation rate constants. For example, for r-i » ri, the probability 
distribution of mRNA levels has a single peak at a low level. As r\ increases, a second peak appears at a 
high expression level with the peak becoming more prominent as r\ approaches r%. At the other extreme of 
parameter values, r± » ra, a single peak at the high expression level is obtained. As r\ changes, the position 
of the peak remains the same, a characteristic of binary response. Stochastic binary gene expression refers to 
a bimodal distribution of mRN A/protein levels and the bifurcation from a unimodal to a bimodal distribution 
occurs in the parameter regime n < 1 and r% < 1. 

In the case of miRNA-mediated regulation of gene expression, the most prominent contribution to binary 
gene expression is obtained when both a\ and 04 + 05 are < in Eq. (14), i.e., when the following inequalities 
are satisfied: 



(18) 



new ^Jnew 

with (fr = — 7 "- 2 tA7, " j) , One notes that the parameter region in which stochastic binary gene expression 
may be obtained is expanded from that (^"- < 1 and ^ < 1) in the case of unregulated gene expression. 
Figure 3 shows the steady state mRNA distribution p(m) versus m when a\ is < and 0,4 — 05 < for the 
parameter values (a) r± < 1, r>i < 1 and (b) r± > 1, r>i < 1. The other parameter values are 7„ elu = 0.1 and 
Inew = 0.2. For each case, the values of A are A = 0.001, A = 1 and A = 20 respectively. We observe that 
increasing (decreasing) values of A disfavours (favours) bimodality. For 7* eji , = 0.2, A = 0, i.e., the case of 
unregulated gene expression, binary gene expression occurs in the parameter regime r\ < 1, r-i < 1 but not in 
the regime r% > 1, r-i < 1. Figure 3 demonstrates the enhanced possibility of binary gene expression in the case 
of miRNA-mediated regulation of gene expression. One notes that if the parameter A is set to zero in Eq. (16), 
the parameter regime in which binary gene expression is observed is expanded to include all values of the ratio 
r\ = — and r% = < 1. The parameter A (= fc °//+7 m * ) is analogous to the dissociation constant for the 
formation of the bound mRNA-miRNA complex with A = signifying infinitely strong binding. We now provide 
a physical understanding of the enhanced occurrence of binary gene expression when a miRNA-regulated gene 
expression parameter, say A, is changed. 

As mentioned in the Introduction, molecular titration can generate ultrasensitive input-output responses. 
The origin of ultrasensitivity lies in the sequestration of an active component in an inactive complex through 
binding to an antagonist [U[TUJ|30]. The simple mechanism involves the kinetic scheme 

X + Y^XY (19) 

where X represents the active molecule, Y its antagonist and XY the inactive complex. Let Kd be the dissociation 
constant associated with the complex XY. The total concentrations of X and Y are given by Xt = [X] + [XY] 
and Yr = [Y] + [XY]. When Yr >> Kd, the X versus Xt plot exhibits ultrasensitivity with the threshold set 
by Yr- For Xt < Yr, almost all the X molecules are part of the inactive complex XY so that the concentration 
of free X molecules is very low. At the threshold point Xt — Yr, the sequestration of X molecules is no longer 
dominant, so that a small increase in Xt can give rise to a large increase in X. The sharpness of the ultrasensitive 
response in the simple example is controlled by the ratio [TU]. In the present study, X and Y represent the 
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Figure 4: (a) Steady state concentration of mRNA to versus total mRNA concentration too (Eq. (20)). The 
parameter Q t h has the value 2 with A = 1CP 5 (dash-dotted line), 0.001 (dashed line), 0.01 (thin solid line) and 
0.1 (thick solid line). Marked ultrasensitive response is obtained for very low values of A. In (b), (c) and (d), 
the p (to) versus to plots are shown for to > Qth (b) , to ~ 9th (c) and m < t h (d) . In the last case, the 
second peak position is at a value very close to zero. 
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mRNA and the miRNA respectively. When 04 + 05 is < in Eq. (14), a dominant singularity occurs at the 
point m = d 5 . with the expression for d$ given in Eq. (16). In terms of the original parameters, d 5 is given by 

d 5 = l[(m - A - e th ) + [(mo - A - 9 th f + 4m A] 1 / 2 ] (20) 

where mo = — , 9th = — = 7m "" T and A is given by the expression in Eq. (6). The magnitude of mo is a 
measure of the total amount of mRNA. When the effective dissociation constant A — > 0, we obtain 



m = d 5 = l - [m - 9 th + ^ [m -9 th ) 2 ] (21) 

= if m < 9 th (22) 
= m if m > 9 t h (23) 

The constant 9 t h is proportional to the total amount, mix, of miRNAs and sets the threshold of an ultrasensitive 
steady state response. Figure 4 (a) shows the steady state m versus mo plots (Eq. (20)) on a semi-logarithmic 
scale for the values of A = 10~ 5 (dash-dotted line), A = 0.001 (dashed line), A = 0.01 (thin solid line) and 
A = 0.1 (thick solid line). One finds that the ultrasensitive behaviour is prominent only when A is close to zero. 
The value of 9 t h is fixed at the value 2. Figures 4 (b)-(d) show the steady state p(m) versus m curves (Eq. 
(14)) for the parameters values 7* eM , = 0.2, 7„ eu , = 0.1 and mo greater than 9 t h (b), almost equal to 9 t h (c) and 
less than 9 t h (d). In the last case, the second peak at m = d 5 (a 4 + a 5 < in Eq. (14)) occurs very close to the 
peak atm = (ai < 0). 

The inequalities in Eq. (18) can be rewritten in a simpler form 

/ ''" I , -^-<1 (24) 



7e//l 7e//2 

with j effl = ^ + 7 m and -f eff2 = 7,„(1 + y^f). 

The rate constants j e ffi and 7 e //2 can be interpreted as effective mRNA decay rate constants by noting the 
following. From Eq. (9), for z = 0, i.e., when the target gene is in the inactive state, 

dm 7*m ... 

-7T = — r - 7m m = / m 

dt m + A 

On linearising / (m) around the steady state value m = 0, one obtains 

dm .7* . 

— = -{— +m)m= -Jeffi m (25) 

Similarly, for z = 1, i.e., when the target gene is in the active state, /(m) = k m — —7m m. Again, linearising 
/ (m) around the steady state value m = dg, one gets 

dm , , 1 — f . , ,_„. 

= k m - 7 m (1 + Y^Tf' m = km ~7eff2 m (26) 

Eqs. (23) and (24) are identical in forms to the equations obtained in the case of unregulated gene expression. 
In the latter case, however, 7 e //i = 7e//2 = 7m- The origin of binary gene expression in the parameter regime 



< 1, < 1 now has a clear physical interpretation. In the inactive state of the gene, the activation 

rate constant has a lower value than that of the effective mRNA degradation rate constant 7 e //i so that the 
accumulated protein level can decay to the steady state level m = before the next transition to the active state 
of the gene occurs. Similarly, in the active state of the gene, the deactivation rate constant has a magnitude 
lower than that of the effective mRNA degradation rate constant Jeff 2- The gene is thus in the active state for a 
sufficiently long time so that the steady state mRNA level m — d^ (Eq. (20)) is attained. The parameter regime 
(Eq. (22)) corresponds to binary response in gene expression as changing the k a , kd values does not alter the 
peak positions but only the amplitudes. When both ^ and le i 2 are > 1, the activation/deactivation rate 
constant is larger than the effective mRNA degradation rate constant, so that a transition from the inactive to 
the active state and vice versa occurs before the mRNA level can attain the value m = in the inactive state, 
and the value m = d 5 in the active state. As a result, the probability distribution p (m) is unimodal, with the 
peak position at an intermediate mRNA level. One also obtains a graded response to changing values of the 
rate constants k a and kd, i.e., the position of the peak shifts in a graded manner. In the steady state expression 
for p (m), for a\ < and + a 5 < 0, the singularities occur at m = and m = d§. The other terms add a 
small background correction to the distribution. Figure 5 shows the plot of p (m) versus m in the case when 
a\ < and — < 0, (a 4 + a 5 > 0). In this case, though a bimodal distribution is obtained for n, r 2 > 1 
and A small, the second peak is associated with quite a broad distribution ((a) and (b)). When the magnitude 
of A is increased from A = 0.002 (b) to A = 0.2 (c), a unimodal distribution is obtained. 
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Figure 5: p (to) versus m plots in the steady state for a% < and — a$ < 0. The values of the parameters 
are quoted in individual boxes. The other parameter values are 7„ et0 = 0.1, 7* letu = 0.2. For low values of A, a 
bimodal distribution is obtained ((a) and (b)) with the second peak more prominent for the lower value of A. 
When A is increased from the value 0.002 in (b) to the value 0.2 in (c), the bimodality in destroyed and one 
gets a unimodal distribution. 



4 Conclusion 

Theoretical modeling studies combined with experiments have been spectacularly successful in uncovering novel 
features of cellular phenomena. One such feature is that of binary gene expression, in which the distribution 
of mRNA/protein levels has two prominent peaks. The possible origins of binary gene expression are several of 
which three principal mechanisms have been the focus of recent studies [20] : positive feedback-based [121 [T3j , 
emergent bistability |15l [2T] and purely stochastic [181 l2"2l 1231 l2~il [231 [231 12"7] . The first two mechanisms create 
the potential for the coexistence of two stable expression states and noise- induced transitions between the states 
give rise to the bimodal nature of the distribution of mRNA/protein levels. In the case of binary gene expression 
with a purely stochastic origin, there is no bistability in the deterministic limit. Experimental observations of 
stochastic binary gene expression have been reported in a number of studies |241 [231 [27] . in agreement with 
theoretical results. The large number of studies carried out so far on binary gene expression consider the 
expression to be unregulated or transcription-factor regulated. The issue of binary gene expression in the case 
of post-transcriptional regulation of gene expression has been mostly confined to the investigation of models in 
which both transcriptional and post-transcriptional (e.g., miRNA-regulated) modes of regulation are considered 
and the transcriptional regulation involves a cooperative positive feedback. In this paper, we have discussed two 
possible scenarios for miRNA-mediated regulation of gene expression, one deterministic and the other stochastic, 
and demonstrated the existence of binary gene expression in each case. In the first case, a nonlinear protein 
decay term, together with a noncooperative positive feedback term generate bistability in specific parameter 
regimes. In the second case, the conditions for obtaining stochastic binary gene expression are obtained and 
one finds that the parameter region in which binary gene expression occurs is more extended in comparison 
with the cases of unregulated and transcription factor-regulated gene expression. In the latter two cases, slow 
transitions between the inactive and active gene expression states are responsible for binary gene expression. 
The conditions — < 1 and < 1 can be reexpressed as > 1 and > 1 with T a , T4 being the average 

lifetimes of the active and inactive gene expression states (T a ~ Td ~ ^-) and T m is the average lifetime 
of the mRNAs. In the case of miRNA-mediated regulation of gene expression, ultrasensitivity generated by 
molecular titration plays a key role in the generation of binary gene expression. The conditions 7 fc ° < 1 and 

^ 4 a < 1 (Eq. (22)) set new time scales in the characterisation of the 'slowness' of the transitions between 
the gene expression states. This results in an expanded parameter regime in which binary gene expression 
is observed. The theoretical predictions made in the present study could be tested in experiments involving 
natural and synthetic gene circuits. 

Binary gene expression generates phenotypic heterogeneity in a cell population with identical genetic makeup 
and exposed to the same initial conditions and environment. Microorganisms adopt a number of strategies for 
coping with stressful situations like environmental fluctuations, nutrient depletion and application of antibiotic 
drugs. One such strategy is the creation of phenotypic heterogeneity so that the whole population does not suffer 
the same fate when subjected to stress [28, 29, 30 . Recent experiments on E. coli [28, 29J, B. Subtilis [30] and 
M. Smegmatis [311 132| have demonstrated the advantages of phenotypic heterogeneity when microorganisms are 
subjected to stress. The regulation of gene expression by miRNAs is known to be activated in a number of cases 
in which a cell population is subjected to stress [TJ[51[3]- Experimental evidence of phenotypic heterogeneity in 
the form of binary gene expression due to regulation by miRNAs would be of significant interest in the context 
of the response of cell populations to stressful conditions. 
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